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ABSTRACT 

We show that spatial correlations in a stochastic large scale velocity field in an 
otherwise smooth intergalactic medium (homogeneous comoving density) su- 
perposed on the general Hubble flow, may cause a 'line-like' structure in QSO 
spectra similar to the population of unsaturated Lya forest lines which usually 
are attributed to individual clouds with 10^^ ^ iVni ^ 5 x 10^^ cm^^. There- 
fore there is no clear observational distinction between a diffuse intergalactic 
medium and discrete intergalactic clouds. It follows that the Hl-density in 
the diffuse intergalactic medium might be substantially underestimated if it 
is determined from the observed intensity distribution near the apparent con- 
tinuum in high resolution spectra of QSOs. Our tentative estimate implies a 
diffuse neutral hydrogen opacity tqp ~ 0.3 at 2: ^ 3 and a current baryon 
density Higm — 0.08, assuming a Hubble constant Hq = 70 km s^^ Mpc~^. 
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The intergalactic medium [IGM] may be probed by observing absorption features in the 
spectra of QSOs. The continuum on the blue side of the HI Lya emission line is expected 
to be depressed as compared to that on the red side, since Lya absorption in the diffuse 
intergalactic medium leads to an apparently continuous absorption due to the general cos- 
mological expansion (Scheuer 1965, Gunn & Peterson 1965). If this effect can be measured 
accurately, it allows in principle to estimate the intergalactic density of neutral hydrogen, 

n{m). 

In their original paper Gunn & Peterson [GP] estimated the depression in the spectrum 
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of 3C9 to be of the order of 40 per cent and n(HI) to be of about 6 x 10~^^ cm~^ a.t z = 2 . 
This low value of n(HI) is usually taken as an indication for a very high degree of ionization 
of the IGM. 

However, the measurement of the true depression caused by the smoothly distributed 
intergalactic hydrogen is hampered by numerous absorption lines blueward of the Lya emis- 
sion line, which make it difficult to locate the continuum level correctly and to remove the 
contribution to the depression from individual absorbers accurately. There have been sug- 
gested a few methods to overcome this problem and to measure the diffuse neutral hydrogen 
absorption (for a review see e.g. Fang & Crotts 1995). The conclusion is that either 'there 
is no evidence for any Gunn-Peterson effect up to the highest redshifts observable in quasar 
spectra' ( Giallongo et al. 1995), or that tgp = 0.115 ± 0.025 at z ~ 3.4 (Fang & Crotts 
1995). In any case the effect is small, indicating a very low density of neutral hydrogen in 
the diffuse intergalactic medium. 

It is worth pointing out that in quantitative estimates the IGM is usually considered as 
a thermal gas and the possible influence of a hydrodynamical velocity field superposed on 
the general Hubble flow is ignored. Commonly the individual Lya lines being interpreted 
as arising in discrete intergalactic clouds at different redshifts are evaluated by means of 
a Voigt profile fitting analysis (see e.g. Lu et al. 1996). In this procedure there are three 
free parameters, the column density A^hi, the Doppler parameter 6, and the redshift Za., to 
be determined from the observational data. Implicitly it is also assumed that any hydrody- 
namical velocities within the individual clouds may be accounted for in the microturbulent 
approximation in which any bulk motions are treated as completely uncorrelated (i.e. the 
velocity correlation length Iq = 0). 

Thus, in the standard analysis it appears observationally possible to clearly distinguish 
between the Lya absorption in the diffuse intergalactic medium leading to a true GP-effect 
and that in the individual clouds leading to the Lya forest lines. As will be shown below, 
this conclusion is misleading and is a consequence of the simplifying assumptions underlying 
the common approach. As a result the diffuse neutral hydrogen opacity mentioned above 
may have been considerably underestimated. 

With respect to the saturated Lya forest hues, corresponding to clouds with 10^^ ^ A^hi ^ 
10^^ cm^^ it has been pointed out recently (Levshakov & Kegel 1996), that the quantitative 
interpretation may be changed substantially when one accounts for a stochastic velocity field 
with finite correlation length within the individual clouds. In particular it has been shown 
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Figure 1. MC simulations of QSO absorption spectra for three random realizations of the stochastic velocity field along a 
given line-of-sight. Panels a, b, and c show a portion of the spectrum blueward of Lya emission for = 3 computed with the 
parameters go = 0.5, Ho = 70 km s~^ Mpc~^, and n(HI) = 4 x 10~^^ cm~^. In addition to the mesoturhulent spectra,, results 
obtained under the condition that the bulk motions have zero correlation coefficient and ctjv^h > 1 (panel a), or have a finite 
correlation coefficient but (Tt/vth ^ 1 (panel b) are shown in order to demonstrate the microturbulent solution and the case of 
purely thermal broadening, respectively. Panel c shows both a high-resolution and a low-resolution spectrum with FWHM = 
0.12 A and 80 A, respectively. The number density of the absorption features in this portion of the simulated QSO spectrum 
is N{z) ~ 500 which is in a good agreement with observations (cf. Kirkman & Tytler 1997). 

by means of Monte-Carlo [MC] simulations (Levshakov, Kegel & Mazets 1997) that such a 
velocity field can lead to a rather structured line profile with several subcomponents, which 
in the traditional analysis would be attributed to individual clouds while in the underlying 
model a homogeneous density was assumed. 

This latter result led us to investigate in which way the GP-effect is modified when a 
stochastic velocity field with a finite correlation length is superposed on the general Hubble 
fiow with a homogeneous comoving density. In fact our results show (Fig. 1) that such 
fields can give rise to a 'line-like' structure of absorption features similar to the observed 
unsaturated Lya forest lines (10^^ ^ A^hi 5 x 10^^ cm~^). 
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Figure 2. An example of MC simulations, (a) - One realization of the stochastic velocity field (in units of the rms velocity 
ct) versus the cosmological redshift Zcosm- (b) and (c) - The corresponding local absorption coefficients for fixed Aa = 4743 A 
and A;, = 4753.5 A, respectively (thin curves, correlation length I ^ 0), and for comparison k\ for a purely thermal gas (thick 
curves, i = ). (d) — A portion of the resulting absorption spectrum. The simulated intensities at the fixed wavelengths Aa and 

Aj, arc marked by dots and the corresponding cosmological redshifts Za and 25 are shown on panel (a). Panel (e) illustrates 
the effect of the contribution from different volume elements to the absorption at a given value of A when the peculiar velocity 
field (panel a) is superposed on the general Hubble flow. 

To describe our mesoturbulent model (/q > 0) the following set of the parameters has 
to be specified : the neutral hydrogen number density n(HI) at cosmic time t, the kinetic 
temperature Tfem, the ratio of the rms value of the line-of-sight stochastic velocity dispersion 
at to the thermal velocity Vth, the correlation length Iq, and the parameter go- We take the 
Hubble constant in units of Hq = 70 km s~^ Mpc~^ and the cosmological constant A = 0. 
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Figure 3. A section of the simulated Lya forest spectrum shown in Fig. la (thin curve) fitted by the superposition of 18 Voigt 
profiles (thick curve). The average deviation of the two curves (xj per degree of freedom) is indicated. 



The stochastic velocity field is assumed to be governed by a stationary Markov process (see 
Levshakov & Kegel 1994, 1997 for details). 

The value of at may be estimated from the observed distribution function of peculiar 
velocities of galaxies. Recent observations indicate, for example, large-scale motions with 
f^gai ^ 700 km s~^ on scales up to 70 Mpc (Raychaudhury & Saslaw 1996). However, if 
the peculiar bulk motions are produced by gravity they should be smaller at high redshifts. 
An order of magnitude oi at at z — 3 may be estimated from the linear theory of the 
gravitational instability which leads to at oc {l + z)~^^^ when the density parameter Q = 2qo 
is equal to 1 (see e.g. Zel'dovich & Sunyaev 1980). For our MC simulations we therefore 
adopt at = 300 km s~^, as well as qo = 0.5 and Zg = 3. The choice of the kinetic temperature 
of Tkin — 10^ K is based on the arguments given by Levshakov & Kegel (1996). For the 
correlation length Iq we chose a value of 24 Mpc. Assuming the velocity fiuctuations to be 
'frozen' into the IGM and just expanding along with the general expansion, we then have 
lz^k/{l + z) =6 Mpc. 

An order of magnitude of n(HI) may be estimated from the average value of Da - the 
continuum depression between HI Lya and Ly/? emission lines in low dispersion spectra 
with FWHM~ 80 — 100 A [Da = (1 — Fobs/Fint), where Fobs and Fint are the observed and 
intrinsic fluxes per unit wavelength in the QSO rest frame]. At 2; = 3 one finds Da — 0.25 
according to the data compiled by Kennefick et al. (1995). 

We will consider in some detail a narrow spectral range of AXobs = 30 A between the Lya 
and P emission lines. For this range we assume that the contribution to Da = 0.25 comes in 
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general from the cumulative effect of weak absorption features, and we will show that in this 
case the Doppler effect caused by the peculiar velocities superposed on the Hubble flow leads 
to the observed picture - instead of a characteristic smooth GP-trough one observes a 'line- 
hke' structure. We emphasize that this does not rule out the possibility of a clumpy IGM, 
but our results indicate that the interpretation of the current high resolution observations 
of QSOs is not unique and therefore the classical analysis may yield misleading parameters 
of the IGM. 

Having specified the free parameters we can simulate a QSO absorption spectrum for 
a random realization of the velocity field. All our synthetic spectra were calculated with a 
signal-to-noise ratio of 100, an instrumental resolution of 8 km s^^, and a pixel size of 2 
km s~^ to match the data obtainable with the HIRES spectrograph of the Keck I telescope. 

We calculate ta, the total optical depth at wavelength A, following GP, the essential 
difference being that in our model the local absorption coefficient kx{z) is a stochastic 
variable (see Fig. 2b, c) which we estimate by means of our MC technique (Levshakov, Kegel 
& Mazets 1997). 

Fig.l shows three mesoturbulent spectra calculated for three different reahzations of 
the velocity field. For comparison and in order to validate our procedure, we calculated in 
addition tx for zero correlation coefficient / = and at/vfh > 1, and for a finite correlation 
coefficient but (Jt/vth ^ 1- For both cases one expects to obtain a classical GP-trough. 
The first case corresponds to the microturbulcnt limit, while in the second case thermal 
broadening dominates. The results of these tests are shown in Fig. la and lb, respectively, 
and marked by 'GPmicro'- 

With n(HI) = 4 X 10~^^ cm~^ at 2; = 3 our microturbulent results are found to fit the 
continuum depression Da ~ 0.25. The corresponding GP-optical depth tgp is about 0.3. 
- Fig. Ic shows a high resolution mesoturbulent spectrum and for comparison the same 
spectrum degraded to a resolution of 80 A. The latter resembles the microturbulent spectra 
in panels a and b. - Commonly one attempts to estimate n(HI) in the diffuse IGM, i.e. 
inbetween Lyo; clouds, by measuring the intensity distribution near the apparent continuum 
level in high resolution spectra of QSOs (see e.g. Fang & Crotts 1995). The corresponding 
regions in our mesoturbulent spectra are labeled by 'GPmeso' in Fig- 1- It is clearly seen that 
instead of tgp — 0.3 one may obtain only an upper limit or a very small value for tgp in 
this case. 

In order to illustrate these calculations. Fig. 2a shows one random realization of the 
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Table 1. 'Lya cloud parameters' derived from the spectrum shown in Fig. 3 by fitting Voigt profiles 



No. 




AT (HI), cm-2 


6, km s ^ 


1 


4743.76 


6.93(11) 


16.10 


2 


4744.26 


6.65(11) 


12.95 


3 


4744.78 


1.84(12) 


16.56 


4 


4745.34 


2.16(12) 


18.50 


5 


4745.82 


2.03(12) 


13.13 


6 


4746.33 


7.21(12) 


18.99 


7 


4746.82 


7.89(12) 


17.32 


8 


4747.43 


1.28(13) 


17.88 


9 


4747.92 


6.60(12) 


22.78 


10 


4749.18 


6 78(^12") 


37.12 


11 


4749.78 


1.53(13) 


30.43 


12 


4750.84 


4.02(1,3) 


24.54 


13 


4751.91 


1.15(13) 


23.01 


14 


4752.50 


1.10(12) 


21.12 


15 


4753.18 


6.33(12) 


19.30 


16 


4753.72 


3.00(12) 


23.51 


17 


4754.50 


2.32(13) 


26.43 


18 


4755.60 


3.63(12) 


17.05 



stochastic velocity field along a given line of sight. - The peculiar velocities have the effect 
that different volume elements may contribute to the absorption at the same value of A. 
This is shown in Fig. 2b, c by the thin curves. For example, the contribution to comes 
from individual volume elements spread over 1 Mpc. For comparison also the absorption by a 
purely thermal gas (no peculiar correlated velocity field) is shown by the thick curves. Finally 
Fig. 2d shows the resulting absorption spectrum and Fig. 2e the effect of compensating the 
general Hubble expansion by hydrodynamic perturbations. 

In order to stress the difference in the interpretation of weak Lya forest lines we fitted a 
part of the simulated spectrum (4743 < A < 4756 A) shown in Fig. la in the usual way by 
considering 18 individual Voigt profiles. The result is shown in Fig. 3 and the corresponding 
parameters are listed in Table 1. As can be seen from the figure the fit is satisfactory, the 
derived 'physical parameters' are in the ranges 6.6 x 10^^ < iV(HI) < 4.0 x 10^^ cm^^, and 
12.95 < b < 37.12 km s~^ (note the thermal width in our model corresponds to 12.89 km 
s~^). So, one can obtain a typical spread of 'physical parameters' for the Lya forest lines 
(cf. e.g. Kirkman & Tytler 1997), while a homogeneous model was assumed. 

Thus, our calculations offer a new interpretation of the observed Lya forest lines and 
explain why high resolution spectra show only a very small GP-cffcct. Our results show that 
at least some of the Lya forest lines may be caused by absorption in the diffuse medium 
between the intervening clouds giving rise to saturated lines. This implies that in the diffuse 
IGM n(HI) may be substantially higher than estimated from the upper limits for tqp in 
previous studies. For example Giallongo et al. (1995) give a value tqp — 0.01 ±0.03 for PKS 

© 0000 RAS, MNRAS 000, 000-000 



8 S.A. Levshakov and W.H. Kegel 

2126-158 {ze = 3.3) which imphes n(HI) ^ 2 x 10~^^ cm~^. This is an order of magnitude 
below the value we assumed in our model. - If one wants to derive from n(HI) the total 
baryon density, one has to know the degree of ionization of the IGM. With the assumptions 
that (i) the ionization of the IGM is in equilibrium, (ii) the ionizing continuum can be 
described by a power law with spectral index = —1.5 and a flux at the hydrogen Lyman 
limit ~ 3 X 10-2^ ergs s"^ cm-^ sr"^ Rz'^ (see e.g. Bechtold 1994), (iii) n(HI) = 4 x 10"^^ 
cm~^, (iv) Tkin = 10^ K, and (v) all the IGM density is in the form of hydrogen, an order of 
magnitude estimate leads to a current baryon density of flicu — 0.08, i.e. it is about 8 per 
cent of the mass needed to close the Universe. 

Taking this at face value, we conclude that Qigm — 0.08 is consistent with the most 
conservative upper limit on baryonic mass ft^ ^ 0.06 from Big Bang Nucleosynthesis (e.g. 
Fields & Ohve 1996). It follows that the density of the IGM may be considerably larger than 
the observed density of luminous matter 0.003 ^ ^^lum ^ 0.007 (e.g. Jcdamzik et al. 1995), 
and the baryon density in the Lja clouds, f^Ly — 0.002 — 0.003 (Lanzetta et al. 1991). 

While our model is highly idealized (smooth density), the results clearly show that the 
formation of the GP-depression is intimately related to the formation of narrow absorption 
lines. Thus, there is no simple way to distinguish observationally between the diffuse IGM 
and intervening clouds as was suggested by the classical interpretation of the GP-effect and 
the Lja forest lines. 
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